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Abstract 

We argue that existing measurements of ep collisions at HERA -in which an ener- 
getic photon is made via a QED 'Compton' subprocess- can provide rather detailed 
information on the photonic parton density of the proton. This function and its 
deviations from Bjorken scaling should be measurable, allowing for an interesting 
test of the theory. The photonic distribution function and its gluonic counterpart 
should show a strikingly different evolution with momentum scale. 
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1 Introduction 



The ep collisions in which there is an energetic photon in the final state can be 
classified according to their distinctive kinematical features. We are concerned in 
this note with the Compton subprocess ej — > ej, where the initial photon is coupled 
to the proton and is (almost) on-shell. Progress at HERA JIJ and the subsequent 
prospects for improved measurements induce us to revisit, update and extend earlier 
theoretical work on this subject [@]-@. 

Figure 1 shows the lowest-order 'Compton' Feynman diagrams, along with our 
notation for the kinematical variables. At the HERA collider, events of this type 
have a particularly clean and distinctive signature, as there are only an electron^] 
and a photon in the final state, with little or even no observable hadronic activity. In 
addition, to lowest order in perturbation theory, the (large) transverse momenta of 
the electron and the photon approximately balance each other. Conversely, applying 
these distinctive features as selection criteria allows one to extract, to the leading 
log approximation, Compton events from the ensemble of radiative corrections [[F]. 




s = (p + lf = x 7 s u = (l- I'f = -Q 2 t={l-k) 2 « -s-u 

Figure 1: Lowest-order Feynman diagrams for Compton scattering in ep collisions. 

We are interested in the ep — > X Compton process at relatively large momen- 
tum scales Q 2 ; we refer to it as Deep Inelastic Compton Scattering (DICS). The 

1 We call both electrons and positrons 'electrons'. 



1 



experimental study of DICS at HERA offers the unique possibility of measuring the 
photon-parton content of the proton and the corresponding longitudinal-momentum 
distribution or 'structure' function j(x, Q 2 ). In the 'parton model' approximation 
in which the initial photon is on-shell and collinear with the proton, the relation 
between the DICS differential cross section and j(x, Q 2 ) is ||: 

d 2 a(s,x,y) f 1 dz d 2 a n ^ n (xs, xj z,y) 
dx dQ 2 

where 



f dz d 2 a e ^ e ^{xs,x z,y) 2 

It d(x/z)d Q 2 7( ^' g) ' (1) 



d 2 a^(s,x,y) 2na 2 em 1 + (1 - yf 
d&dQ* = ^ I^y— 5{1 - X) - (2) 



A measurement of j(x, Q 2 ) would acquire particular interest when compared 
with that of its non-Abelian counterpart: the gluon distribution function g(x,Q 2 ). 
The predictable Q 2 evolution of these functions ought to be very different, their 
difference representing a very clean test of the non-Abelian nature of the gluon. 
More specifically, the self-coupling of gluons makes g(x, Q 2 ) evolve with Q 2 , at small 
x, much faster than 7(2, Q 2 ). The evolution of 7(2, Q 2 ) with Q 2 is by itself an 
interesting test of the QCD-based parton picture, which we argue can be performed 
at HERA with the existing or soon to be gathered statistics. 

One can in principle, and sometimes in practice, distinguish two types of contri- 
butions to DICS: 'pseudo-elastic' ep — > e^p, and 'inelastic' ep — > e-yX with X ^ p. 
From a parton model point of view, in which the struck photon in Fig. 1 is viewed 
as a proton constituent, the distinction between pseudo-elastic and inelastic sub- 
processes is very artificial. This is best understood by comparing conventional deep 
inelastic scattering with the Compton process, as we do in Fig. 2. Deep inelastic 
scattering, Fig. 2a, results in two final-state jets: the 'current' or struck quark jet 
and the 'spectator' jet of target fragments. The invariant masses of these jets (to 
the extent that they can be ascertained without ambiguity) are small; the unavoid- 
able colour reconnection between the struck quark and the target fragments is in 
general a soft process. At high Q 2 a successful recombination leading to an elastic 
event ep — > ep is a rare ocurrence. Similarly, the photon and the target fragments 
in the Compton process of Fig. 2b would even more rarely recombine into a proton, 
to result in a literally elastic scattering event ep — > ep. There is nothing very spe- 
cial about the 'pseudo-elastic' ep — > e^p channel, which ought to be quite common, 
since the invariant mass distribution of the target fragments will, as in deep inelastic 
scattering, tend to be low. In DICS the total invariant mass of the proton's frag- 
ments (the hadrons as well as the outgoing photon) will, again as in deep inelastic 
scattering, be large in general. 

The pseudo-elastic contribution to DICS can be easily worked out in terms of 
measured proton form factors ||. In the current state of our understanding of 
QCD, it is simply impossible to predict the contribution of final hadronic states of 
invariant mass Wh > m p . To study the feasibility of measuring 7(2, Q 2 ) and its Q 2 
evolution in DICS we must, as other authors have done before us, proceed 
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(a) 



(b) 



Figure 2: Diagrams for deep inelastic conventional and 'Compton' scattering. 

to make a series of conjectures. Mercifully, these conjectures are not relevant to the 
Q 2 evolution, which is a solid prediction of the standard model. 



2 Pseudo-elastic DICS 



The contribution of the pseudo-elastic channel ep — > epy to 7(37, Q 2 ) can be explicitly 
written down in terms of the proton's elastic form factors ||. The usual 'electric' 
and 'magnetic' form factors are empirically well fit as dipoles: 



G E (t) 



Define the quantities: 



1 



[1 - 1/(0.71 GeV 



2\12 



, G M {t)~ 2.79 G E (t). 



(3) 



G 2 (t) - (t/4m 2 )G 2 Af (t) 
1 — t / 4m^ 

to express the pseudo-elastic contribution to 7(2, Q 2 ) as: 



(4) 



7dW 



a, 



-x 



2tt ./ tl t 



d ±l2 



H x {t)+H 2 {t) 



(5) 



where the limits on the virtual photon mass t are: 



(s + m p )(s [l — x]+ m p ) ± (s — m p )J(s [1 — x] + m 2 ,) 2 — 4sm 2 



(6) 
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The tempting analogy between photons and gluons as partons in the proton 
breaks down at various points. One of them is that colour confinement and conser- 
vation preclude the existence of a strict coloured analogue to a quasi-elastic channel: 
there are no QCD elastic form factors. A neutral and spinless hadron, such as a K^, 
would have a very small quasi-elastic contribution to its photonic structure func- 
tion Jk(%, Q 2 ), satisfying in this respect a QED/QCD analogy more closely than 
a proton does. But the proton carries a long-range photon field and j e i( x ) is an 
important contribution to its 7(2, Q 2 ). Another departure from the photon/gluon 
analogy is that, in DICS, it is justified to work to leading order in a em , and to this 
order the quantity 7 e ;(a;) of Eq. (|) is independent of Q 2 . 

At an ep collider, a separation of elastic and inelastic Compton events is difficult 
but possible. A fraction of elastic events could be tagged by the very forward 
'diffraction' detectors. Inelastic events with proton fragments at sufficiently large 
Pt could be caught by larger-angle detectors. Although the elastic contribution 
can be ascertained with accuracy, we see no particular interest in singling it out 
experimentally except, perhaps, as a calibration and/or luminosity check [I]]. 



3 Inelastic DICS 



In principle one could build-up the complete function 7(2, Q 2 ) by adding to the 
quasi-elastic contribution all resonant |J and non-resonant final hadronic states, 
and their interferences. Alternatively one could directly guess, one way or another, 
an inclusive or 'continuous' non-elastic part of 7(2, Q 2 ). If this guess is based on a 
parton picture wherein the photon 'constituent' is emitted by one of the quarks in 
the proton |3|, [|, |(J , the addition of the continuous and resonant contributions may 
be double-counting, as it would certainly be in e + e~ annihilation || and arguably 
be in deep-inelastic scattering || . In these latter cases the data provide the decisive 
proof that there is a 'duality' between continuous and resonant contributions: adding 
them is double-counting. In DICS, the information on 7(2, Q 2 ) is at the moment 
too sparse to help decide on this interesting question. 



Having to make a guess, we choose to make one that is correct in that aspect 
of j(x, Q 2 ) for which the QCD prediction is unequivocal: the Q 2 evolution. Two 
diagrams for the contribution of quark-gluon scattering to high-p^ photon produc- 
tion are shown in Fig. 3. Notice that, with the substitution of gluons for photons 
and of a parton-quark for an electron, they are identical to the Compton diagrams 
of Fig. 1. This implies that DICS is 'factorizable' in the same sense as 'Drell-Yan' 



scattering or the production of high-p^ photons or jets in hadronic collisions | ID| . In 
turn, this means that 7(2, Q 2 ) satisfies a 'QED evolution' equation [11]. To lowest 
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order in a em and a s 0, and in an obvious notationf]: 

di(x, Q 2 ) a em f 1 dy 



dlnQ 2 2ir 



I - E < p m (-) kv, Q 2 ) + v(y, Q 2 )} , (7) 
j x y w j 



where the quark-to-gauge-boson splitting function is: 

P M(Z) = l±l±-ll . (8) 

A function 7(2;, Q 2 ) satisfying Eq. (|7|) is automatically correct to the leading loga- 
rithmic order in QED. This is also the order to which it is permissible to single out 
Compton scattering from the other radiative processes in ep collisions. 

We have chosen the lepton-to-lepton squared momentum transfer Q 2 as the scale 
in Eq. (0). A scale p^(e, 7) or either of the quantities s, —t defined in Fig. 1 would 
be equally reasonable. To leading order of perturbation theory there is no way to 
favour any particular choice. 




The non-Abelian generalization of Eq. (0) is the gluon-evolution equation ||12|| : 



dg(x,Q 2 ) = a s (Q 2 
dlnQ 2 



2n 



'dy 

y 



E ^ p m ( - ) [i(y, Q 2 )+i(y, Q 2 )\ +p 99 [-) g(v, Q z 



(9) 

where the (lowest-order) gluon-to-gluon splitting function P gg is given in []T3]. An 
interesting challenge to experiment is to observe the different Q 2 evolutions of 7 and 
g, as predicted in Eqs.(|7p and (Bp. 

To illustrate the experimental feasibility of a measurement of 7(2;, Q 2 ) and its 
evolution, we must guess that function explicitly. For the guess to be consistent with 
QCD, it must be compatible with Eq. (|7|). We follow |J in writing 7 = 'jei + lin and 
figuring out 7i„(x,Q 2 ) by evolving an input "fi n (x, Qq) with the use of Eq. (0). We 
also follow H in assuming that 7j n vanishes at a low Qq = 0.5 GeV; it builds up at 

2 The 0(a em ) modifications of the quark and gluon evolutions induced by a non- vanishing 
7(21, Q 2 ), affect the evolution of the latter only at 0{a 2 em ). 
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Q 2 > Qq via evolution. This is admittedly as arbitrary for a photon constituency 
as it would be for its gluon counterpart; the proton being a bound state of charged 
coloured objects, there is no reason for it to be 'made' of only quarks at any scale, 
even at leading twist and leading order of perturbation theory. We differ from || in 
our use of a more recent set [|n| of parton densities^. 



We deal with momentum scales comparable to the proton mass and we should be 



making target mass-corrections distinguishing Bjorken's a;-scaling |TJ| from Nacht- 
mann's ^-scaling [15]. As it turns out, 7(2, Q 2 ) is only measurable at values of x 
small enough for this distinction not to be important. 



4 Experimental details and expectations 



In [|TJ , a first measurement of DICS at HERA was presented; at that time the main 
aim was the use of Compton events as an independent luminosity monitor. The 
roughly 400 reported events were insufficient to study differential distributions and 
to provide a measurement of 7. Moreover, for many of these events it was not pos- 
sible to tell which of the electromagnetic showers was the electron and which was 
the photon. Since then, about two orders of magnitude more luminosity have been 
accumulated. In addition, silicon trackers have been installed in the backward re- 
gion, and charge identification is now possible with much higher efficiency We show 
that, given these improvements, a lot of information on 7(2, Q 2 ) can be extracted 
from existing data. 

To match the experimental situation, we use the following parameters and cuts, 
the latter akin to the ones used for event selection in the quoted HI analysis [[J: 



• The beam energies are E e = 27.5 GeV and E p = 820 GeV. We slightly over- 
estimate event rates by using the current value of £ = 36.5/pb for the total 
luminosity^ collected by HI. 

• The electron and the photon are both seen by the central detector, i.e. they 
have 0.05 rad < 6 eri < it — 0.05 rad (as usual, angles are measured with respect 
to the proton beam direction). 

• At least one of the electromagnetic clusters has an energy of more than 8 GeV, 
and both of them have at least 5 GeV. The latter requirement guarantees a 
high detection efficiency. 

• Other cuts in are automatically satisfied at the level of our leading-order 
calculation. The total visible electromagnetic energy is always larger than 18 

3 This leads to an increase of the conjectured j(x, /1 2 ) at small 2, relative to the result in [||. 
4 Our numbers can be easily rescaled to the luminosity collected after the installation of the 
silicon trackers. 



6 



GeV (the final e + 7 energy equals E e + x y E p > 27.5 GeV). Except in the 
forward region there is no additional (hadronic) cluster with more than 2 GeV 
energy. The acoplanarity angle A0 = n — \<p e — 7 | is below 45°. 

We also demand that no momentum scale be small: for the quantities defined in 
Fig. 1, —t, s, Q 2 > 1 GeV 2 ; for the (equal and opposite) transverse momenta of the 
final-state electron and photon, pff > 1 GeV. 

The first relevant question concerns the number of DICS events, with the quoted 
integrated HERA luminosity. In Fig. 4 we show, in bins of log 10 x, the event numbers 
that survive the above selection criteria. Clearly, they should suffice to measure this 
x-distribution, a point previously emphasized in 0-0. The electron and photon 
detection efficiency (and presumably the ability to distinguish them) is high within 
the applied cuts. We shall thus estimate the statistical errors as the square root of 
the event number per bin. To measure x, e/7 identification is unnecessary, since 
x = x 7 may be determined from the requirement that the sum of the energies of the 
two electromagnetic clusters be equal to the combination E e + x^E p . Contrawise, 
for the determination of Q 2 , an e/7 distinction is indispensable. 
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Figure 4: Event rates for the DICS at HERA. The cuts applied are as described in the 
text. The lowest-x bin contains all events with log 10 (x) < —4.25. 



In Fig. 4 the dashed (dash-dotted) histogram depicts the number of DICS events 
found by integration of Eq. (0) using a scale s (— t ) instead of Q 2 . The scale 
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dependence is a substantial theoretical uncertainty, but is not devastatingly large. 
The fact that the elastic part is scale-independent helps make the total theoretical 
expectation fairly stable. The dotted histogram in Fig. 4 shows the event rates 
predicted on the basis of quasi-elastic scattering only: 7 = 7^. It should be an easy 
task for experiment to confirm, or infirm, the existence of the substantial inelastic 
contribution that we foretell. 

To illustrate the experimental extraction of 7(2, Q 2 ) we translate the information 
in Fig. 4 into a statement on the accuracy of the measurement. To this end, we 
evaluate j(x, Q 2 ) at the statistical averages (x),(Q 2 ) determined from the event 
sample used in Fig. 4. We assume that in each bin the error in 7 is only statistical. 
The result for /a em as a function of log 10 (x) is shown in Fig. 5. For comparison, we 



also show the contribution of j e i. In practice a possible method |16| to translate the 
data into a measurement of 7 is based on the bin-by-bin evaluation and subsequent 
iteration of the expression: 

i"~K*>, (Q 2 )] - # * CT M q^] • VK*>, « 2 >] • (I") 

where #et>., MC[7 toy ] is the expected (theoretical or MonteCarlo) number of events, 
based on a trial function 7 toy . If r y toy is not too far off the true structure function 7, 
the iteration of Eq. ( |10|) converges fast. A reasonable guess is 7* oy = 7 e ;. To see how 
well this choice works, we make use of it on the right-hand side of Eq. (|T0D, along 
with our results of Fig. 4 as the 'data sample'. The x r j meas / a em so determined is 
also shown in Fig. 5. Even this 'zeroth' order iteration is accurate. 

To study the detectability of the Q 2 dependence of 7, we take the sample used for 
Figs. 4 and 5 and bin it additionally in Q 2 . The result is shown in Fig. 6. There is 
an increase of xj(x, Q 2 )/ot em with Q 2 that should be observable in most x-bins. We 
also display the Q 2 dependence of the leading-order GRV |L3| gluon density, which 
satisfies Eq. (|9|). To facilitate the comparison of the relative variations of 7 and 
g, we renormalize xg(x,Q 2 ) in each subplot so as to coincide with xj(x, Q 2 )/a em 
at the lowest accessible Q 2 . The gluonic function g evolves much more strongly at 
small x than its photon counterpart 7. Only at x > 10 -2 does the Q 2 -evolution of 
7 overtake. The statistics appear to be insufficient to explore in detail the region 
x > 0.07, where g starts to decrease with Q 2 , while, for 7, Eq. (|7]) predicts an 
increase for all x. 



5 Conclusions 



We have shown that it is possible to extract j(x, Q 2 ) -the function describing the 
photon constituency of protons- from existing data. Even its Q 2 -dependence should 
be already observable, allowing for an interesting test of a combination of QED 
and QCD. The photon and gluon distribution functions should show a strikingly 
different Q 2 evolution. 
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Figure 5: Expected statistical accuracy of the determination of X7((x), (Q 2 ))- The 
numbers indicate the average (Q 2 ) (in GeV 2 ) for each x-b'm. The dotted line shows 
the purely elastic spectrum, while the dashed one corresponds to >y meas 'extracted' via 
Eq. (pH|) and as explained in the text. 

We know so little about the photonic structure function 7 that the possible im- 
provements on its theoretical understanding would be premature, since the extra 
precision they would bring about is surely at a more refined level than the current 
uncertainty. The corrections to next-to-leading order in a s would be the most inter- 
esting, as they should help control the dependence of the predictions on the choice 
of variable representing the momentum scale. As measurements of DICS materialize 
-and we hope they soon will- these improvements should become timely. 
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Figure 6: Expected statistical accuracy of a measurement of the Q 2 dependence of 
x"f(x, Q 2 )/a em in various a>bins. The dashed line displays the Q 2 dependence of the 
LO GRV [[13[] gluon density, which has been normalized in each plot so as to coincide 
with xj(x, Q 2 )/a em at the lowest accessible Q 2 value. 
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